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Abstract. We have used a Monte Carlo radiative transfer code to produce edge-on images of dusty galactic 
disks, allowing a fraction of the dust to be distributed in clumps. Synthetic images of edge-on galaxies have been 
constructed for different amount s of dust, distributions of clumps and fractions of dust in clumps, following the 
formalism of Bianchi et al. (2000| ). We have also considered models with stellar emission embedded in the clumps. 
The synthetic images have been fitted with ana lytical models made with smooth distributions of dust, adopting 
the procedure developed by Xilouris et al. (1999) to fit optical images of real edge-on galaxies. We have compared 
the parameters determined by the fit with the input parameters of the models. For the clumping distributions 
adopted in this paper, the neglect of clumping results in underestimating the amount of dust in a galaxy. However, 
the underestimation is never larger than 40%. 
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^ ' 1. Introduction 



Optical extinction by interstellar grains complicates the 
study of the stellar and dust content of a spiral galaxy. 
For moderate to high optical depths, information on the 
intrinsic properties of each of these components is en- 
tangled in the observed surface brightness distribution of 
the galaxy. The usual way of deriving these properties 
is by comparing the observations with radiative transfer 
calculations in a dusty galactic model. A certain degree 
of complexity is needed in the description of the geo- 
metric distributions of a galactic model. Radiative trans- 
fer models with simple geometries have been proven to 



provide equivocal results (Disney, Davies, & Phillipps 
198£). ^Furtbermnre^ scattering muHt hp taken into a.c 



and thus can not be easily included in fitting procedures. 
For such procedures an approximate treatment of the ra- 
diative transfer is more appropriate. 

The first approach of such a procedure was introduced 
by Kylafis & Bahcall (1987, hereafter KB), where verti- 
cal profiles of a model galaxy were fitted to observations 
of NGC 891. In KB's work, scattering is calculated up 
to the first order and an approximation is introduced for 
the higher order contributions while the distribution of 
the dust is assumed to be smooth and exponential along 
the radial and vertical directions. A similar approximation 
for the treatment of scattering is adopted by Silva et al. 
(1998). Other works prefer an exact solution for scatter- 
ing within simple geometries (Bruzual, Magris, & Calvet 



count, given tho high albodn nbsprvpd for Oalartir dust 1988; [Corradi, Bcckman, fc Simonneau 1996| ; |Xu fc Buat 



( iGordfin Cal^ntti W\tt IQQTl) 



1995) 



Using the KB method, the radiative transfer equation 



A radiative transfer code using the Monte Carlo (MC) 
technique can, in principle, accommodate any stellar and can be solved for a wide variety of geometries (Byun et al 
dust d stribution and provide an accurate calculation of 1994). Xilouris et al. (1997;1998) improved the original 



the surface brightness for a given model (Witt, Thronson 
& Capuano 1992; Bianchi, Ferrara, & Giovanardi 199t ; |de 
Jong 1996| ). However, the MC method is time consuming 
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idea of KB and implemented a technique to fit the sur- 
face brightness distribution of edge-on spirals. In edge-on 
galaxies, the effects of extinction are maximized and it 
is possible to separate the stellar and dust components. 
The method has been successfully applied to a sample 
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of seven edge-on galaxies (Xilouris et al. 1999) conclud- 
ing that late- type spiral galaxies have a moderate opacity, 
with a mean face-on optical depth in the B-band tb ~ 0.8. 
As in most of the models used for the description of spiral 
galaxies, Kilouris et al. (199£) adopt smooth exponential 
distributions to describe the stellar and dust disks. 

Real galaxies exhibit a wide variety of inhomogeneities, 
like spiral arms, bars, clumps. While it is possible to in- 
clude these structures in the solution, when fitting an ob- 
served image it is desirable to deal with the simplest pos- 
sible description, in order to limit the number of model 
parameters. Complex models can then be used to test 
the reliability of the description obtained with the sim- 



ple models. This is done, for example, in Misiriotis et al. 
(2000|)] where synthetic edge-on images of model galax- 
ies with a spiral structure are fitted with a smooth ex- 
ponential model. They find that plain exponential distri- 
butions provide a good description of the galactic disks. 
The derivation of the stellar and dust parameters is only 
slightly affected by the spiral pattern. 

In this paper we perform a similar exercise to study 
the influence of dust clumping on the fit of edge-on galax- 
ies. In most clumpy medium exhibits higher 
transparency than a smooth one of the same dust mass. 
Therefore, the comparison of real images with smooth 
models may result in an underestimation of the dust con- 
tent in a galaxy. Here we use the MC method presented 
by Bianchi et al. (2000, hereafter BFDA) to create im- 
ages of clumpy galactic models. For the parametrisation of 
BFDA, the maximum difference between clumpy and ho- 
mogeneous models can be seen in the edge-on case. Monte 
Carlo models of highly inclined galax ies had been pro- 
duced also by Kuchinski et al. (1998| ), with a different 
parametrisation for the clump distribution. With respect 
to their result, the influence of clumping in the edge-on 
case is larger for the model adopted here (see the discus- 
sion in BFDA). The M C images of BFDA wi ll be fitted 
with the KB model as in Misiriotis et al. (2000 ), aiming at 
determining the efficiency of clumping in hiding the dust 
when a galaxy is seen edge-on. 

Sect. ^ describes the method adopted for the com- 
parison and the MC models we have used to produce 
the synthetic images. A model with exponentially dis- 
tributed clumps developed for this work is presented in 
the Appendix. The results are presented in Sect. || and a 
summary is given in Sect. ^. 

2. Our Method 



Following the procedure described by Misiriotis et al 
[(2000| ), we fit images made with the KB method to the 
MC images of a clumpy model galaxy. Both the MC and 
the KB m odels presented here are com puted for the V- 
band. The Henyey fc Greenstein (1941 ) phase function is 
used with the scattering paramet ers in the V band al bedo 
oj = 0.6 and asymmetry g = 0.6 ( Gordon et al. 1997 ). 

The parameters defining the smooth models are the 
central face-on surface brightness I^,, the radial scalelength 



a* and the vertical scalelength f3i, of the stellar disk, the 
V-band central face-on optical depth ry, the radial scale- 
length ad and the vertical scalelength /3d of the dust disk. 
The total intrinsic luminosity L-i, is equal to 47r/^aJ while 
the dust mass Md is propotional to Tya^. The results from 
the fits are compared to the input parameters of the MC 
images and are analyzed as a function of the fraction of 
dust locked into clumps fc- Additionally, we study the 
importance of scattering in the fitted models. 

The disk parameters adopted for the MC models are 
those used in BFDA. For the stellar disk, we have adopted 
a ratio between radial and vertical scalelengths a^//3* — 
11.5. The smooth dust disk has the same radial scalelength 
as the stellar one and a vertical scalelength /3d = 0.4/3^,. 
In both the simulation and the fitting procedure the dust 
distribution has been truncated at 4.6 radial scalelengths 
along the radial direction and at 4 vertical scalelengths 
along the vertical direction. 

In the clumpy images, Ty defines the total mass of 
dust in the model, i.e. for the same ry the clumpy and 
the smooth models have the same dust mass. In partially 
clumpy models, a fraction fc of the total dust mass is 
distributed in clumps, while the remaining is distributed 
in a smooth disk of central face-on optical depth tv(1 — /c). 
The number of the clumps Nc is determined as described 
in BFDA. 

To cover the full range of ry and fc we considered the 
following optical depths : rv=0.5, 1, 2 and 5. For each 
optical depth we produced images with /c= (smooth 
case), 0.25, 0.5, 0.75 and 1.0, resulting in a set of 20 im- 
ages. Several different sets can be produced depending on 
how we distribute the dust clumps. We have decided to 
produce 3 sets. 

In the set named "molecular" , the dust clumps follow 
the same distribution as the molecular gas in the Galaxy 
and their characteristics are based on the Giant Molecular 
Clouds. Details on these models can be found in BFDA. 
In this set, the /c=l case is not considered because it 
would mean that all the dust is distributed in a molecular 
ring and an exponential model is not reccommended to 
reproduce such a distribution. 

A second set of models named "exponential" was pro- 
duced with the dust clumps following the same exponen- 
tial distribution as the smooth dust disk while the stars 
were smoothly distributed. The exponential models are 
described in detail in the Appendix. 

In the third set of models named "embedded" the 
clumps are distributed exponentially, while a fraction of 
starlight /(,,„b=0.5 was embedded in the dust clumps. For 
the case of /c=0, there can be no /emb=0.5 because when 
there are no clumps there is nowhere to embed the em- 
bedded starlight. 

The MC images of the clumpy model galaxies were 
constructed to be as similar as possible to observations. 
The images cover 151 pixels along the major axis from the 
galactic center to the border of the image, while the radial 
scalelength of the stellar and the dust disk is equal to 50 
pixels. For each image, photons are gathered in a finite 
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Fig. 1. Contour plot of a synthetic edge- 
on image from a Monte Carlo smooth 
model with tv=1. Overplotted are 



isopho tes from a Kylafls 
(1987|) type model. 
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inclination band (Bianchi ct al. 1996). which was chosen 
to be 1°. The number of photons for the MC simulations 
was chose n so that a realist ic signal to noise ratio was 
achieved ( Bianchi ct al. 1996 ). 

The fitting was performed as in [Misiriotis et al. (2000| ), 
by minimizing the sum of the squares of the differences be- 
tween the MC and the KB images. In Fig. |l| we present a 
contour plot of an edge-on image from a MC smooth model 
(/c=0) together with the isophotes of the same model im- 
age, produced using the KB method. The coincidence of 
the MC and KB model images demonstrates the consis- 
tency of the two methods. 



3. Results of the fits 

3.1. Testing of the codes and the effects of scattering 

Since the two radiative transfer methods (BFDA and KB) 
used in this work have been developed independently and 
with different algorithms, we checked for their consistency 
in several ways. We first produced images without scatter- 
ing and without clumps using the MC method and fitted 
them with the KB method. The fits were able to retrieve 
all the parameters with a scatter smaller than 4%. 

We then analyzed the smooth models (/c = 0) de- 
scribed in the previous section. Two fits were produced 
for each MC image. In the first, the fitting model is com- 
puted without scattering, while in the second scattering 
is taken into account. This way we can quantify both the 
importance of scattering in fits of edge-on galaxies and the 
approximation of KB in their treatment of scattering. 

In Fig. H we show the parameters J* and ry whose be- 
havior changes the most when passing from an absorption- 
only fit to a fit including scattering. It is well-known 
that scattering reduces the effective absorption of a dusty 
medium, thus resulting in an increased surface brightness. 
When scattering is not included (dotted lines in Fig. 
the fitting procedure compensates for this by reducing the 
optical depth of the dust disk and by increasing the central 
surface brightness of the stellar disk. The underestimate 
of the dust content is the dominant effect, because the 
stellar luminosity is significantly constrained by the parts 
of the galaxy that lie outside of the dust lane, where the 
absorption is negligible. Absorption-only fits result in an 
underestimate of the face-on optical depth ranging from 
10% for Ty=0.5 to 25% for Ty=5. As for the central face- 



on surface brightness (/*), the maximum discrepancy is 
an overestimation by 10% for tv=0.5 

When scattering is included, the fitted parameters are 
very close to the input values (solid lines in Fig. The ry 
plot confirms the va lidity of the approximat ion for scat- 
tering introduced by Kylafis fc Bahcall (1987|) , at least for 
low-to-moderate optical depths. Even for Ty=5, the de- 
rived optical depth is only 15% larger than the input value. 
All the other parameters derived from fits with scattering 
were found to vary within less than 4% of the input values. 

3.2. The effect of clumping 

The quality of the fits is poorer when synthetic images 
from clumpy models are analyzed. The sum of the squares 
of the differences between input data and fitted models in- 
creases by about an order of magnitude. Furthermore, in 
many cases it was possible to find a few equally acceptable 
fits, corresponding to local minima. While some parame- 
ters of KB models may change from minimum to mini- 
mum, we found that the derived dust mass and intrinsic 
luminosity are nearly the same. 

In Fig. ^ we have plotted the derived normalized dust 
mass and intrinsic luminosity, as functions of and of ry , 
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Fig. 2. The derived and ry normalized to the input val- 
ues as functions of the optical depth of the model. Results 
for fits with and without scattering are shown. 
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Fig. 3. Normalized dust mass (left column) and normalized intrinsic luminosity (right column) derived by fitting the 
Monte Carlo simulations of Section as functions of the fraction /c of dust locked up in clumps. The values have 
been normalized to the input dust mass and intrinsic luminosity. Results are presented for models with optical depths 
Tv= 0.5 (dotted lines), 1 (dashed), 2 (dot-dashed) and 5 (solid). The plots in the top row show the results for the 
clumpy model of BFDA with fraction of embedded stellar emission /omb=0; the middle row refers to the models with 
clumps distributed exponentially and /emb^O; the bottom row to exponentially distributed clumps, but with /cmb=0.5. 
Results for the smooth models (/c = 0) are plotted in both the top and middle row . 



for the "molecular" models, (BFDA; top row); for the "ex- 
ponential" models (middle row) and for the "embedded" 
models (bottom row). 

The derived dust masses, normalized to the input val- 
ues, are shown in the left column of Fig. The trend 
of the dust mass with /c does not depend on the type 
of clumping or on the optical depth of the model, the 
amount of detected dust decreases with increasing /c. For 



the two sets of models without embedded stellar emission 
( "molecular" and "exponential" ) , the dust mass is under- 
estimated by less than 10% for /c=0.25 to about 20-30% 
for /c=0.75. For the "exponential" set for /c=l, the un- 
derestimate reaches values of 40%. For the "embedded" 
models the fitting procedure recovers the input dust quite 
well. All the dust is detected for /c=0.25 and only 20% is 
missed for /c=0.75. 
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The results of the fitting procedure for the intrinsic, 
dust-free, luminosity L^, arc shown in the right column of 
Fig. ^ again normalized to the input value. As for the dust 
mass, the trends with /c for different models are very sim- 
ilar and do not depend sensitively on the value of ry. The 
fitting model overestimates the luminosity with increasing 
/c. For fc—0.75 and /omb = 0, the derived luminosity is 
about 20% larger than the input value. 

The overestimate of the luminosity drives the fitting 
procedure to increase the dust mass in order to obscure 
the additional starlight. If the luminosity was not over- 
estimated then even less dust would be needed to fit the 
synthetic images. A simple major-axis analysis, i.e. a com- 
parison of the radial profiles along the dust lane suggests 
that the smooth model detects only half of the dust mass 
(BFDA) and leads to the conclusion that clumping causes 
an underestimate of the face-on optical depth. 

This is not the case for the fits shown here. The fitting 
procedure tries to interpret the reduced extinction caused 
by clumping by decreasing the dust radial scalelength aa 
and keeping the optical depth high (actually, to a value 
similar to the input one), rather than decreasing Ty and 
keeping ad fixed. The reason for this may be that the 
surface brightness in models with clumping departs from 
that of a smooth exponential disk, when all the image, 
and not only a single profile, is fitted. Given the general 
trend of the retrieved data, it does not seem probable that 
we have systematically hit local minima of the fit. 

We have shown that the trend of dust mass underes- 
timation is nearly independent of the optical depth of the 
model. This may suggest that the same percentage of dust 
is missed when fitting a galaxy at different wavelengths. 
Therefore, clumping should not affect the determination of 
the extinction law from fits of a galaxy in different bands. 
As long as the dust properties do not vary significantly 
with wavelength, as in the optical range, it is possible to 



of an edge-on clumpy galactic model is analyzed through 
a simpler (and more controllable) smooth model. As ex- 
pected, using a smooth model when fitting images of edge- 
on galaxies leads to an underestimation of the dust content 
of the galaxy. For the distribution of clumps adopted in 
this paper, about 20-30% of the dust mass is missed by 
the fitting, the value depending more on the fraction of 
dust located in clumps rather than on the details of the 
clumps distribution or the total dust mass of the model. In 
addition, clumping alters the surface brightness distribu- 
tion with respect to that of a smooth exponential model, 
thus making the fit more difficult. This may explain the 
overestimation of the intrinsic luminosity in models with 
clumping. 

We stress again here that our work is limited to the 
edge-on case, for which convergent fits can be obtained 
with the KB technique. The results may also depend on 
the clump distribution and on the geometrical parameters 
adopted here. However, it is interesting to note that other 
works find a small influence of clumping in edge-on galax- 
ies. Kuchinski et al. (1998) produced radiative transfer 



models of disk galaxies including clumping. The distribu- 
tion of clumps is different from the one adopted here: the 
clumpy medium is characterised by a costant filling fac- 
tor throughout the whole disk and by a density contrast 
with the homogeneous medium ( Witt fc Gordon 1996| , see 
also BFDA for a comparison with the present models). 
Comparing colour gradients across the extinction lane of 
highly inclined galaxies with the model results, they find 
that the inferred opacity of the dust disk is largely insen- 
sitive to the difference between clumpy and homogeneous 
dust distributions. 

Recent analyses of dust emission at A > 100//m sug- 
gest that the amount of dust in spiral galaxies is larger 
than wh at derived from fits of edge-on surface brightness 
profiles (Bianchi et al. 2000 ; Popescu et al. 2000 ; Misiriotis 



interpr et the absence of a trend with ry as an absence of 
a trend with A. 



4. Summary 

We have used a Monte Carlo code for the radiative trans- 
fer in a clumpy galactic model to produce synthetic edge- 
on images of spiral galaxies in the V-band. The images 



et al. 2001 ). Clumping has been imputed to as a possible 



cause for this discrepancy: dust in star-forming regions 
may contribute dominantly to the Far Infrared emission 
without being seen through its extinction effects. However, 
the underestimation produced by the clumpy distribu- 
tions adopted in this paper is always smaller than 40%, 
unable to explain the magnitude of the effect. Different 



have t len been analyzed using the fitting procedure of 



Xilouris et al. (1999| ), that uses the method of Kylafi; 
fc Bahcall (1987| ) to compute the radiative transfer in 
smooth galactic models. We first checked for the consis- 
tency of the two independent codes, by analyzing smooth 
MC images. Despite the differences in the solution of the 
radiative transfer, it was possible to retrieve the param- 
eters of the MC models within 4% of their original value 
and we confirmed the importance of scattering in the ra- 
diative transfer simulations. The dust mass is significantly 
underestimated when scattering is not included in the fit- 
ting routine. 

We have then studied the influence of clumping on 
the derivation of the parameters, when a synthetic image 



dust/clump distributions may be needed. Popescu et al 
(2000D and [Misiriotis et al. (200 iD find that the long wave- 



length emission observed in two galaxies (namely, NGC 
891 and NGC 5907) can be explained by including a sec- 
ond dust disk, at least as massive as the dust disk causing 
the observed extinction lane and thin enough to escape 
detection through the technique of surface brightness fit- 
ting. 
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Appendix A: Models with exponentially 
distributed clumps 

In BFDA it was assumed that all the material associated 
with the molecular phase in a spiral galaxy is in Giant 
Molecular Clouds. Radiative transfer models were pro- 
duced assuming that dust clumps are distributed as the 
molecular gas in the Galaxy, i.e. in a ring-like structure. 
However, most of luminous, face-on, late-type galaxies, 
and a good fraction of early-type spirals, show a ra dial ex- 
ponential distrib ution for the H2 column density ( Young 
& Scoville 1991 ). For the sake of completeness, we briefiy 
present here a set of models with an exponential distribu- 
tion for the dust clumps. The clumps distribution has the 




K [a.] 

Fig. A.l. Major-axis edge-on profiles for t\j~2 mod- 
els with molecular (top panel) and exponential (bottom 
panel) clumping, for /c=0.25,0.5 and 0.75 (from bottom 
to top, solid lines). As a reference, we also show the profile 
for a TV =2 smooth model (dotted lines). 

same scalelengths as the smooth dust disk. All the other 
parameters are the same as in BFDA. 

The effects of clumping in these models are not very 
different from those obtained adopting the Galactic H2 
distribution. Exponentially distributed clumps make the 
models less transparent, but only slightly. For example, in 
the Tv=2 model with fraction of gas in clumps fc=0.75, 
15% of the radiation is absorbed, while in BFDA it was 
12%. We recall that in a smooth model 19% of the ra- 
diation is absorbed. For models with embedded stellar 
emission, the difference is even smaller, because of the in- 
creased contribution of the absorption within each clump. 
The number of clumps and their opacity are independent 
of the assumed clump distribution and are the same as in 
BFDA. 

As in BFDA, the largest differences between homoge- 
neous and clumpy simulations are for the edge-on profiles. 
In Fig. AT we show the edge-on major-axis profiles for the 
BFDA models and for the models with exponential clump- 
ing, for rv=2, /c=0.25,0.5,0.75 and fomb=0. The profile 
for a smooth rv=2 disk is shown as a reference. Because 
of the exponential distribution, the profiles of the new 
clumpy models are similar to those for smooth exponen- 
tial disks, independently of /c and fomb- Instead, BFDA 
models show the effects of the ring-like clumping distribu- 
tion, even when fcmb=0. 



